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Cold and ultracold polar molecules with nonzero electronic angular momentum are of great interest 
for studies in quantum chemistry and control, investigations of novel quantum systems, and precision 
measurement. However, in mixed electric and magnetic fields, these molecules are generically subject 
to a large set of avoided crossings among their Zeeman sublevels; in magnetic traps, these crossings 
lead to distorted potentials and trap loss from electric bias fields. We have characterized these 
crossings in OH by microwave-transferring trapped OH molecules from the upper |/; M = + §) 
parity state to the lower |e; +§) state and observing their trap dynamics under an applied electric 
bias field. Our observations are very well described by a simple Landau-Zener model, yielding insight 
to the rich spectra and dynamics of polar radicals in mixed external fields. 

PACS numbers: 37.10.Pq,33.80.Be,33.20.-t 



Cold polar molecule experiments have made a remark- 
able transformation over the past decade, progressing all 
the way from the first demonstrations of basic motional 
control of molecules [1-3] to studies in ultracold chem- 
istry [4], polar collisions [5-7], and precision measure- 
ment [8, 9[. Closed-shell molecules at low temperatures 
and long ranges present nearly featureless electric dipoles 
[10, 11], yielding universal scattering behavior and the 
prediction of generic dipolar condensates [12] or crystals 
[13]. At shorter ranges, the quantum statistics and chem- 
ical nature of the specific species comes to the forefront 
[4, 14]. Polar radicals behave in even richer fashions: 
topological crystals [15], half-integer vortices generated 
by conical intersections [16], and spin-dependent chem- 
istry [17] have been predicted, and the utility of using 
magnetic fields to trap radicals while performing electric- 
dipole-dependent studies has already been demonstrated 
[14, 18, 19]. Motivated by these predictions, several 
groups are pursuing the production of ultracold radicals 
such as RbSr [19, 20] and LiYb [21]. 

While the combination of electric and magnetic dipoles 
enables fascinating new phenomena, it also comes with 
substantial complications. In particular, simple linear 
Zeeman spectra are converted to tangles of interwoven 
avoided crossings by the application of remarkably small 
transverse electric fields. In this Letter, we report the 
experimental observation of these crossings in the polar 
radical OH and their dramatic impact on dynamics in 
a magnetic trap. In its 2 n 3 / 2 ground electronic state, 
stationary OH has a magnetic dipole moment of 2^b 
(where jig is the Bohr magneton) and an electric moment 
of 1.67 D (0.65 a. u.); rotational averaging reduces each 
of these in the laboratory frame to fi[ e ^^ — l J -ijfj + \), 
where J is the total angular momentum, M is its lab- 
fixed projection, and £l is the magnitude of J's projec- 
tion on the internuclear axis. Each rotational level of 
OH contains two Zeeman manifolds of opposite parity, 



|e; M) and |/; M); with no applied fields, these states 
are split by a A-doublet coupling of 1.667 GHz. In an 
applied electric field, |e) and |/) are not strictly good 
quantum numbers. However, they are still useful as labels 
of a state's Stark character: |e) states are strong-electric- 
field-seeking while |/) states are weak- field-seeking. 

A representative Hamiltonian for A-doublet molecules 
in mixed electric and magnetic fields has been derived 
in [22]. While realistic systems must be solved nu- 
merically, the basic structure of the combined H to tai — 
H zero— field + H stark + H zeeman Hamiltonian is simple: 
the zero-field part generates diagonal matrix elements in 
a suitable basis (e. g. Hund's case A fine-structure plus 
parity basis |jMf2e) for OH, where e = e or / is the 
J-relative parity label), as does the Zeeman part if one 
takes the space-fixed Z axis to lie along the magnetic 
field B. The off-diagonal elements are then purely from 
the Stark portion. The component of the electric field 
E parallel to B generates matrix elements which pre- 
serve M and change the parity e, while the transverse 
components change M by ±1 as well as e. As can be 
seen in the Zeeman spectrum of OH (Fig. 1(a)), all 
states in the lower parity manifold except for the ab- 
solute ground state have Zeeman-induced crossings with 
one or more of the upper parity states: in an E-field, 
these real crossings are coupled by the off-diagonal Stark 
elements and become avoided. Working in the pertur- 
bative limit, the crossings Xm of Fig. 1 have widths 
S+i/2 ot \E sin 0eb\, <5-i/2 oc \E 3 sin 2 9eb cos0eb\, an d 
<5-3/2 \E 3 sin 3 6 eb\, where Oeb is the angle between 
E and B. Since the off-diagonal elements only couple 
states of opposite parity, the total perturbation theory 
order must be odd: crossings with total AM even such 
as X_i/ 2 require one extra AM = coupling order so 
that the total coupling changes e as it must. This mix of 
AM = ±1 and AM = matrix elements in turn leads to 
a (sin 9 eb) AM cos 9 eb angular dependence in the cou- 
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Figure 1. (color online) The Zeeman structure of OH with bias 
electric fields: (a) V/cm, (b) 500 V/cm at a 90° magnetic- 
to-electric field angle 6 E b, (c) 5000 V/cm at 65°, and (d) 5000 
V/cm at 90°. The crossings Xm of the je; +§) state with the 
three]/; M — {— f , — §, +§}) states are labeled in (a). 



pling strengths for these even crossings, while the an- 
gular dependence for crossings with AM odd is simply 
(sin0EB) AM ■ Perturbation theory begins to fail, how- 
ever, at relatively small E'-fields (much smaller than the 
polarizing E'-field at zero £?-field) because of the strength 
of the electric dipole coupling and the fact that the Zee- 
man crossings are gapless at zero i?-field. Therefore, for 
all calculations described in this Letter we diagonalize 
the full 8x8 ground state fine-structure Hamiltonian of 
OH [23, 24]. 

We observe these avoided crossings by their dynam- 
ical effects on a sample of magnetically trapped OH 
molecules. Our Stark deceleration and magnetic trap- 
ping system has been described elsewhere [14]. Briefly, 
OH molecules are created in an electric discharge through 
water vapor in a supersonic expansion and then Stark 
decelerated and brought to rest in a permanent mag- 
netic quadrupole trap (Fig. 2(a)) at a temperature of 
~ 70 mK, in their |/; +§) state. By applying chirped 
microwave fields to the trap magnet surfaces, we can 
then transfer the trapped molecules to the |e; +§) state, 
wherein we can probe the three crossings labeled in Fig. 
1(a); application of a DC voltage to the magnets results 
in a bias E- field with the distribution shown in Fig. 2(b). 
The spatial loci of the three crossings in the actual field 
gradient of our trap are shown in Fig. 2(c). The multiple- 
well structure is a result of our use of permanent mag- 
nets: the hole in the center of the magnets allows flux to 
exit both inside and outside the magnet ring, yielding a 
toroidal field minimum between the faces of the magnet, 




Figure 2. (color online, to scale) (a) The magnetic field distri- 
bution in the permanent magnet quadrupole trap. (b). The 
electric field distribution created by applying a potential dif- 
ference to the trap magnets, (c) The spatial location of the 
avoided crossings Xm- (d) The relative angle Oeb between 
electric and magnetic fields in the trap. White denotes 0° 
while black denotes 180°, implying that the avoided crossings 
are widest in the gray regions. 



as well as the central quadrupole trap. The distribution 
of Oeb angles generated by charging the magnets is plot- 
ted in Fig. 2(d). 

As our trapped molecules are moving at velocities < 10 
m/s, only a very small energy gap is needed for each 
molecule to follow the adiabats rather than the diabats of 
the potential surfaces. As can be seen in Fig. 1, the adia- 
bats of all the levels below the "doubly-stretched" |/; +|) 
state are substantially less trapped than the e; +|) di- 
abat - or in some cases, are outright antitrapped. Thus, 
if molecules are transferred to the |e; +§) diabat in zero 
field and then a bias field is applied, some fraction of 
the molecules will leak out of the trap as they follow 
the adiabats rather than Landau-Zener hopping [25] to 
follow the trap diabat. At low fields, the Landau-Zener 
hopping probability is still substantial, and so it takes 
multiple trap oscillations for a molecule to leak out. At 
higher fields, however, effectively all molecules follow the 
adiabats and so the loss rate is simply limited by the 
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where v tra p is the 



trap dynamical time Td yn 
trap oscillation frequency. 

In order to transfer population from !/;+§) to 
|e; +|), we use an Adiabatic Rapid Passage (ARP) tech- 
nique [26] via microwave fields applied directly to the 
trap magnets. As the magnets also experience ±14 kV 



potentials as part of the trap loading sequence, the mi- 
crowave system is AC-coupled to the trap by a pair of 
in- vacuum high voltage capacitors [23]. This capacitive 
coupling stage reduces the ±14 kV bias to a few volts; 
two 1.2 GHz high-pass filters further suppress both the 
near-DC bias and the switching transients. By chirp- 
ing our microwave frequency through the entire range of 
differential Zeeman shifts [8] in our trap over a period 
t~arp <C Tdym we are able to achieve transfer of ~ 75 % 
of the observable trap population. (Since the transfer 
is a 7r transition, only the fraction of Erf \\ B is useful. 
Thus, molecules in the gray regions of Fig. 2(d) see much 
smaller effective intensities than molecules in the white 
or black regions.) 

In addition to causing dynamical trap loss, the avoided 
crossings Xm also reduce the efficiency of our ARP 
pulses: in the vicinity of the avoided crossings, not only 
are the states shifted away from the nominal 1.667 GHz 
transition frequency, but also the two parity states are 
mixed by the electric field. The net result of this mixing 
is that molecules transferred near said crossings may not 
in fact be transferred to the e; +|) diabat but rather 
may end up on the |/; M) diabat and thus be invisible 
to our state-selective detection. Figure 3 displays the 
|e; +|) population as a function of the DC bias field 
applied during the ARP. The bias field leads to a rapid 
drop in transferred population as the |/; +§) — > e; +§) 
transition is Stark-shifted out of resonance; the bias field 
required for this shift increases with the frequency width 
of the ARP pulse. We have performed simulations using 
a numerical solution of the Optical Bloch Equations to 
calculate a position dependent ARP transfer probability 
with zero electric field. We then calculate the position- 
dependent Stark-Zeeman shift by diagonalizing the OH 
Hamiltonian and convolve the resulting frequency distri- 
bution with the ARP bandwidth to determine the total 
fraction transferred. The simulation results shown in Fig. 
3 give a good match to the population measurements, 
showing that the ARP process is well-understood [23]. 
The simulations suggest the presence of a stray E-field in 
the trap region on the order of 100-150 V/cm. 

Applied E-fields also have a profound effect on trap 
dynamics. A loss of ~ 40% of the e-state population over 
Tdyn is observed immediately after the ARP pulse (Fig. 
3 inset), almost certainly due to prompt loss of molecules 
following an adiabat through the A +1 / 2 crossing (which 
is fully avoided in just the observed stray E-field) and 
exiting the trap. At larger bias fields, loss rates up to 
^ T dyn are observed with critical fields again on the order 
of 200-300 V/cm. The observed loss rates (Fig. 4(a)) are 
well-fit by a probabilistic model 
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Figure 3. (color online) Post-ARP e-state population versus 
applied bias field, with 10 (green triangles), 20 (blue squares), 
and 40 MHz (red diamond) ARP chirp widths. Solid lines 
are Monte-Carlo simulations of the population transfer using 
realistic bias and RF field distributions from finite-element 
calculations, a semiclassical treatment of the Stark detunings, 
and a homogeneous stray field as a free parameter, (inset) 
Time-of-flight of e-state population immediately after ARP, 
with no applied bias field. The rapid loss through the X +1 / 2 
crossing is clearly visible: the solid line is a fit to a single- 
exponential loss at a rate of 2r^ n . 



where T 



vacuum 



1 LZ 



is the background vacuum loss rate and 
2^(5 2 M (E eff )) \ 



exp 



MA W (M))(f) 



is the Landau-Zener probability that the molecule stays 
on the e; + §) diabat through the crossing Xm- An- 
gle brackets denote population averages: (^jr) = 
VB ■ v^j is the rate of change of the magnetic field and 

E 2 

stray 



(5 M (E ef f)) = (<5 M (Je 2 + E 2 tray j J is the width of the 
crossing Xm as a function of applied electric field E and 

/ 3/2.3/ 2 _ 3/ 2 . M \ 
^3/2-5/2 3/2-5/2 J 



stray field E stray , while A/i e// (M) - - 

2/ib is the difference between the effective magnetic mo- 
ments of the /; +|) and \e; M) states. The model 
contains only 3 free parameters: T vacuum , E stray , and 
£ « 0.12, which incorporates the errors in our estimates 
of the population-averaged parameters. (While we can 
estimate (^f ) reasonably well from Monte-Carlo trap 
dynamics simulations, (5m {Eeff)) has a very fast vari- 
ation with (0eb) and so is hard to accurately calculate. 
It is therefore responsible for most of £'s deviation from 
unity.) T vacuum can be independently observed from the 
/-state lifetime, while E stra y can be derived from Fig. 3; 
both of these independent measurements are consistent 
with the fit. Td yn is directly measurable from the trap os- 
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Figure 4. (color online) (a) Loss rate of e; M = +§) 
molecules versus applied electric field. The solid red line is a 
fit to a simple Landau-Zener model; the green shaded stripe 
denotes the range of observed loss rates for ]/; +|) molecules. 
The bias field is turned on 5 ms after the ARP pulse is com- 
pleted, and thus well after the prompt loss seen in the inset of 
Fig. 3 is completed, (inset) Widths 5m of the avoided cross- 
ings between e; +§) and /; M) versus applied electric field, 
at their respectively maximizing values of 6eb- <5+i/2 is linear 
with field, while 8_ 1 / 2 and <5_3/2 are cubic, (b) Time-of-flight 
trace of trap oscillations, observed by loading the trap with 
a non-zero center-of-mass velocity. The solid red fit yields a 
trap frequency v trav = 675 Hz. 



cillation frequency v tra p observed in Fig. 4(b). As both 
1/2 an d <5_3/2 are third-order in E, the most striking 
feature of the Landau-Zener model is that it explains the 
extremely steep rise in T: T ~ 1 — exp (— ai? 6 ) where a 
represents the Landau-Zener parameters. 

The success of this model prompts two conclusions. 
The first is that magnetic trapping of the lower parity 
doublet in any molecule is fraught with difficulties: in 
OH, even just the fields from stray charge accumulating 
on the trap mounts is enough to cause unity loss through 
the -X+i/2 channel. (We are currently implementing a 
new trap mount design to minimize the patch charges 
in our future work.) In general, if one wishes to explore 
polar radicals in mixed fields in any states except the 
doubly-stretched ones, either optical trapping or an Ioffc- 
Pritchard trap with a sufficiently large offset B-field as to 
be above all the crossings is necessary; this implies that 
polar radicals in their lower parity doublet are extremely 
challenging candidates for slowing with a magnetic coil- 
gun [27, 28]. More broadly, in precision measurement 
applications it is imperative to take a close look [29] at 
the exact Stark-Zeeman spectrum of the system of inter- 



est, as the spectral distortions from the avoided crossings 
can persist over very large ranges in magnetic field. 

This work was funded by DOE, the AFOSR MURI 
on Cold Molecules, NSF, and NIST. M. Hummon is an 
National Research Council postdoctoral fellow. 
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Derivation of the Stark-Zeeman Hamiltonian 



In order to evaluate molecule-field interaction matrix elements of the form Hi = —ft ■ T (where P represents either 
B or E), it is necessary to transform between the molecule-fixed frame in which the appropriate dipole moment jl is 
defined and the laboratory frame in which T is defined. Following Lara et al. [1], we use the Wigner rotation matrix 
elements V]^ (to) to write 

H i = -a.p=-Y i {-if (^2 Kk h ^-r 

While OH is not fully a Hund's case (a) molecule, it is nonetheless close enough that both the electric and magnetic 
dipole moments may be considered to lie along the internuclear axis and so Hi simplifies to 

H l = -Y,(-^ q Ko (w) Wfc=o.F- g . 

If B and E are not parallel, it is further necessary to rotate the laboratory fields to a combined frame. Without loss 
of generality, we choose Z along B and transform E using the reduced rotation matrix elements d]^ (Oeb) to write 

^ = -VH (co) 4% \B\-J2Ko M &d l * q0 (6 EB ) \E\ . 

<? 

Using Lara et al. 's results for the | JMQe) parity basis yields the matrix elements 
(JMne\Hi\j'M'n'e') = V2 J + ly/2 J' + 1 ^ J J ^ (-1) 

( / . , . J i j' i no \ 



2 ){-M M' jMsS(A + ffe S) 

( —M \ M' ) d -*i" ( 0eb) ^ eE 



+ ( l-ee'(-l) J+J ' +2fr 



where g e is the electron Lande g- factor and /i e is the molecular electric dipole moment. 
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B. Full OH Hamiltonian 

Evaluating the matrix elements described previously for the OH | J = | , M, fi = | , e) ground state manifold yields 
the 8x8 matrix used for combined Stark-Zeeman calculations in this paper: 

/ -#-|f"B B § Me Ecosf) EB -^ fle E ainO EB \ 

- \n B B -^fcEsioSEB l Mc Eco s S EB -§ Me E s i„e EB 

_A + 2 (jbB -§MeE 8 E B - 1 Mc E cos £) E B - ^ Me E sin £> EB 

_^ + | MbB -xfe £ ™»ES -|n e Eco B £) EB 

|l"eS™»EB -^n e Esine EB A - g n B B 

-^(i e E.in8j; B ln e Ecos£> EB - § n e E sin £> E B A - | n B B 

-| Me E Bin£) EB -i|i t E coae EB - ^ Mc E sin E B A + | M B B 

\ -^p e E s in» EJi -g Mc Ecos£) EB A+g MB B y 



C. Microwave coupling circuit design 
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Figure SI: A schematic of the microwave system use to transfer OH between A-doublet levels. 



In order to couple microwave power on to our trap magnets while protecting the microwave system from the ±14 kV 
applied to the magnets during the trap loading sequence, we use a pair of custom in-vacuum capacitors to block the 
DC high voltage plus two more high-pass filters (Mini-Circuits VHF-1200+[2]) to block the AC transients generated 
by the magnet switching. 

Tunable microwave power is generated by an HP 8644A synthesizer, which is programmed to form the ARP pulse 
by voltages supplied to its frequency and pulse modulation inputs. A Mini-Circuits ZFDC-10-5-S+ directional coupler 
allows monitoring of the microwave signal before it is amplified by a GTC GRF5016A power amplifier. The microwaves 
are then coupled through the protection filters and into the vacuum system via an SMA feedthru. Typical power 
levels are 30.4-35.4 dBm (1-3 W) at the feedthru. As the ARP chirp time tarp must be much shorter than Td yn , 
T arp is generally between 20 and 80 /us. 
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D. Numerical simulations of ARP 
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Figure S2: ARP spectra at several different bias J5-fields. Points are data while solid lines are simulation results. The Stark- 
Zeeman broadening due to the avoided crossings is clearly visible. 

As described in the main text, we have performed numerical simulations of the ARP population dynamics. The 
simulation code works as follows: 

1. For each molecule in a Monte Carlo distribution, the Optical Bloch Equations are solved using the local mag- 
nitude of the microwave field (the microwave field is assumed to follow the DC field distribution generated by 
charging the magnet surfaces) and its angle against the local magnetic field. Both field distributions were cal- 
culated using a finite-element package. The final probability of the molecule transferring to the lower A-doublet 
state is then recorded. 

2. For each molecule, the local static E'-field is calculated and then combined with the local B-field to find the 
total Stark-Zeeman detuning. This yields a transfer probability density as a function of microwave frequency. 

3. The probability density function is then integrated over the range of frequencies the ARP chirp sweeps through 
to get a final transfer probability. 

Figure S2 shows several ARP spectra at different bias £?-fields. The agreement between data (points) and simulation 
(solid lines) is good; in particular, the simulations capture the Stark-Zeeman broadening visible at larger fields. In 
a homogeneous field and the absence of avoided crossings, the application of an .E-field would merely push the two 
A-doublet states apart without broadening the transition. The field in the trap region is not fully homogeneous, 
which contributes some broadening, while the curvature of the |e; AI — + §) state near its avoided crossings further 
increases the Stark-Zeeman broadening. In order to fit the simulations to the data, a background stray _E-field on the 
order of 120 V/cm was required, consistent with the field required for the Landau-Zener model. 
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